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Abstract 

The continuous-time quantum walk on the underlying graphs of association schemes 
have been studied, via the algebraic combinatorics structures of association schemes, 
namely semi-simple modules of their Bose-Mesner and (reference state dependent) Ter- 
williger algebras. By choosing the (walk) starting site as a reference state, the Terwilliger 
algebra connected with this choice turns the graph into the metric space with a dis- 
tance function, hence stratifies the graph into a (d+1) disjoint union of strata (associate 
classes), where the amplitudes of observing the continuous-time quantum walk on all 
sites belonging to a given stratum are the same. Using the similarity of all vertices of 
underlying graph of an association scheme, it is shown that the transition probabilities 
between the vertices depend only on the distance between the vertices (kind of relations 
or association classes). Hence for a continuous-time quantum walk over a graph associ- 
ated with a given scheme with diameter d, we have exactly {d + 1) different transition 
probabilities (i.e., the number of strata or number of distinct eigenvalues of adjacency 
matrix) . 

In graphs of association schemes with known spectrum, namely with relevant Bose- 
Mesner algebras of known eigenvalues and eigenstates, the transition amplitudes and 
average probabilities are given in terms of dual eigenvalues of association schemes. As 
most of association schemes arise from finite groups, hence the continuous-time walk 
on generic group association schemes with real and complex representations have been 
studied in great details, where the transition amplitudes are given in terms of characters 
of groups. Further investigated examples are the walk on graphs of association schemes 
of symmetric Sn, Dihedral -D2m and cyclic groups. 

Also, following Ref.[l], the spectral distributions connected to the highest irreducible 
representations of Terwilliger algebras of some rather important graphs, namely distance 
regular graphs, have been presented. Then using spectral distribution, the amplitudes 
of continuous-time quantum walk on strongly regular graphs such as cycle graph C„ 
and Johnson, and strongly regular graphs such as Petersen graphs and normal subgroup 
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graphs have been evaluated. Likewise, using the method of spectral distribution, we 
have evaluated the amplitudes of continuous-time quantum walk on symmetric prod- 
uct of trivial association schemes such as Hamming graphs, where their amplitudes are 
proportional to the product of amplitudes of constituent sub-graphs, and walk does not 
generate any entanglement between constituent sub-graphs. 

Keywords: Continuous-time quantum walk, Association scheme, Bose-Mesner 
algebra, Terwilliger algebra. Spectral distribution. Distance regular graph. 
PACs Index: 03.65.Ud 
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Random walks on graphs are the bases of a number of classical algorithms. Examples include 
2-SAT (satisfiability for certain types of Boolean formulas), graph connectivity, and finding 
satisfying assignments for Boolean formulas. It is this success of random walks that motivated 
the study of their quantum analogs in order to explore whether they might extend the set 
of quantum algorithms. This has led to a number of studies. Quantum walks on the line 
were examined by Nayak and Vishwanath [2], and on the cycle by Aharonov et al. [3]. The 
latter has also considered a number of properties of quantum walks on general graphs[l]. Two 
distinct types of quantum walks have been identified: for the continuous-time quantum walk a 
time independent Hamiltonian governs a continuous evolution of a single particle in a Hilbert 
space spanned by the vertices of a graph [1, 3, 4, 5] , while the discrete-time quantum walk 
requires a quantum coin as an additional degree of freedom in order to allow for a discrete-time 
unitary evolution in the space of the nodes of a graph. The connection between both types of 
quantum walks is not clear up to now, but in both cases different topologies of the underlying 
graph have been studied (see, for example, [2, 6, 7, 8, 9, 10]). 

Different behavior of the quantum walk as compared to the classical random walk have been 
reported under various circumstances. For instance, a very promising feature of a quantum 
walk on a hypercube, namely an exponentially faster hitting time as compared to a classical 
random walk, has been presently found (numerically) by Yamasaki et al. [11] and (analytically) 
by Kempe [12]. Indeed, first quantum algorithms based on quantum walks which offer an 
(exponential) speedup over their optimal classical counterpart have been reported in Refs.[13, 
14]. 

On the other hand, the theory of association schemes has its origin in the design of statistical 
experiments. The motivation came from the investigation of special kinds of partitions of the 
cartesian square of a set for the construction of partially balanced block designs. In this 
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context association schemes were introduced by R. C. Bose and K. R. Nair. Although the 
concept of an association scheme was introduced by Bose and Nair, the term itself was first 
coined by R. C. Bose and T. Shimamoto in [15]. In 1973, through the work of P. Delsarte [16] 
certain association schemes were shown to play a central role in the study of error correcting 
codes. This connection of association schemes to algebraic codes, strongly regular graphs, 
distance regular graphs, design theory etc., further intensified their study. Association schemes 
have since then become the fundamental, perhaps the most important objects in algebraic 
combinatorics. To this regard association schemes have for some time been studied by various 
people under such names as centralizer algebras, coherent configurations, Schur rings etc. 
Correspondingly, there are many different approaches to the study of association schemes. 

A further step in the study of association schemes was their algebraization. This formula- 
tion was done by R. C. Bose and D. M. Mesner who introduced to each association scheme a 
matrix algebra generated by the adjacency matrices of the association scheme. This matrix al- 
gebra came to be known as the adjacency algebra of the association scheme or the Bose-Mesner 
algebra, after the names of the people who introduced them. The other formulation was done by 
P. Terwilliger, known as the Terwilliger algebra. This algebra is a finite-dimensional, semisim- 
ple and is non-commutative in general. The Terwilliger algebra has been used to study P- and 
Q-polynomial schemes [17], group schemes [18, 19], and Doob schemes [20]. 

Here in this paper, we study continuous-time quantum walk on the underlying graphs 
arising from association schemes, by using their algebraic combinatorics structures, namely 
semi-simple modules of their Bose-Mesner and (reference state dependent) Terwilliger algebras. 
By choosing the (walk)starting site as a reference state, the TerwiUiger algebra connected with 
this choice turns the graph into the metric space with a distance function, and hence stratifies, 
the graph into a (d-|-l) disjoint union of strata (associate classes), where the amphtudes of 
observing the continuous-time quantum walk on all sites belonging to a given stratum are the 
same. Since all vertices of underlying graph of an association scheme are similar or they have 
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a constant measure of similarity, therefore the transition probabihties between the vertices 
depend only on the distance between the vertices (kind of relations or association classes). 
Hence for a continuous-time quantum walk over a graph associated with a given scheme with 
diameter d, we have exactly {d + l) different transition probabilities (i.e., the number of strata 
or number of distinct eigenvalues of adjacency matrix). 

In graphs of association schemes with known spectrum, namely with relevant Bose-Mesner 
algebras of known eigenvalues and eigenstates, the transition amplitudes and average proba- 
bilities are given in terms of dual eigenvalues of association schemes. As most of association 
schemes arise from finite groups, hence we have studied in great details continuous-time walk 
on generic group association schemes with real and complex representations, where the tran- 
sition amplitudes are given in terms of characters of groups. Further more, as examples, we 
have investigated walk on graphs of association schemes of symmetric Sn, Dihedral and 
cyclic groups. 

Also following Ref.[l], we have presented the spectral distribution connected to the highest 
irreducible representation of Terwilliger algebras of some rather important graphs, namely dis- 
tance regular ones (since the Hilbert space of walk consists of irreducible module of TerwiUiger 
algebra with maximal dimension). Then using the spectral distribution, we have evaluated 
the amplitudes of continuous-time quantum walk distance regular graphs such as cycle graph 
Cn and Johnson, and strongly regular graphs such as Petersen graphs and normal subgroup 
graphs. Likewise, using the method of spectral distribution, we have evaluated the amplitudes 
of continuous-time quantum walk on symmetric product of trivial association schemes such 
as Hamming graphs, where their amplitudes are proportional to the product of amplitudes 
of constituent sub-graphs, and walk does not generate any entanglement between constituent 
sub- graphs. 

The organization of this paper is as follows. In section 2, we give a brief outline of associ- 
ation schemes, Bose-Mesner and Terwilliger algebras and stratification. Section 3 is devoted 
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to studying continuous-time quantum walk on graphs with known spectrum. In section 4, we 
review continuous-time quantum walk on group association scheme. In section 5, following 
Ref.[l, 22], we investigate continuous-time quantum walk on distance regular graphs via spec- 
tral distribution /j, of the adjacency matrix A. In section 6, we calculate the amplitudes for 
continuous-time quantum walk on some graphs by using the prescription of sections 3, 4, 5. 
The paper is ended with a brief conclusion and three appendices, where the first appendix con- 
sists of studying the method of symmetrization of non-symmetric group schemes, the second 
appendix contains the proof of the lemma regarding the equality of the amplitudes associated 
with the vertices belonging to the same stratum , and the third appendix contains the list 
of some of the finite distance regular graphs with their corresponding spectral distributions, 
respectively. 



2 Association scheme, Bose-Mesner algebra, Terwilliger 
algebra and its modules 

In this section we give a brief outline of some of the main features of association scheme, such 
as adjacency matrices, Bose-Mesner algebra and Terwilliger algebra. At the end by choosing 
the (walk) starting site as a reference state we stratify the underlying graphs of association 
schemes via the relevant Terwilliger algebra connected with this choice. 

2.1 Association schemes 

First we recall the definition of association schemes. The reader is referred to Ref.[23], for 
further information on association schemes. 

Definition 2.1. (Association schemes). Let V he a set of vertices, and let Ri{i = 0, 1, d) be 
nonempty relations on V (i.e., subset oiV xV). Let the following conditions (1), (2), (3) and 
(4) be satisfied. Then the pair Y — {V, {Ri}o<i<d) consisting of a set V and a set of relations 
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{Ri}o<i<d is called an association scheme. 

(1) {Ri}o<i<d is a partition oiV xV 

(2) Ro = {{a, a) -.aeV} 

(3) Ri = R\ for < i < where R\ = {(/?, a) : (a, /?) e Ri) 

(4) Given (q;,/3) G Rk, Pij —\{'yEV: {a,P) e Ri and (7,/9) G -Rj} |, where the constants 
pfj are called the intersection numbers, depend only on i,j and k and not on the choice of 



Then the number n of the vertices V is called the order of the association scheme and Ri 
is called a relation or associate class. 

Let r = (V, R) denote a finite, undirected, connected graph, with vertex set V, edge set R, 
path-length distance function 9, and diameter d ■.—max{d{a,/3) : a,P e V}. For all a,P eV 
and all integer i, we set 



so that Ti{a) = {j3 & V : d{a, j3) — i}. Then T is distance regular graph if Fq, Fi, F^ form 
an association scheme on V. In any connected graph, if /3 e Fj(Q;), then 



{a,p) e Rk. 



Fi = {V, Ri) = {{a, P):a,l3eV: d{a, /3) = z}, 



(2-1) 



ri(/3) C Fi_i(Q;) U Fiia) U Fi+iia). 



(2-2) 



Hence in a distance regular graph, = (for i 0, j is not {i — + 1}) and we set 



cIz^pI, bi^Pl-1,1, Ci^pl^^i, 



(2-3) 



see Fig.l. 



2.2 The Bose-Mesner algebra 

Let C denote the field of complex numbers. By Matv{C) we mean the C-algebra consisting of 
all matrices whose entries are in C and whose rows and columns are indexed by V. For each 
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integer i {0 <i < d), let Ai denote the matrix in MatviC) with (a, /3)-entry 

r 1 if {a, (3) e Ri, 
{A^)a,fS = {a,(3eV). (2-4) 

[ otherwise 

The matrices Ai arc called the adjacency matrices of the association scheme. We then have 
Aq — /(by (2) above) and 

A^, = EpfA (2-5) 

fe=0 

(by (4) above), so Ao, Ai, A^ form a basis for a commutative algebra A of Matv{C), where 
A is known as the Bose-Mesner algebra of y = (V, {-Ri}o<i<d)- 

Then, by using Eq.(2-3) and Bose-Mesner algebra (2-5), for adjacency matrices of distance 
regular graph F, we have 

AiAi = Cj_iA-i + (oi -h- Ci)Ai 

A^Ad = Cd-iAd-i + (ai - bd)Ad, (2-6) 

where A^ is a polynomial in Ai of degree i and A\ is a linear combination of /, Ai, A^. 

Since the matrices Ai commute, they can be diagonalized simultaneously (see Marcus and 
Mine [24]), that is, there exist a matrix S such that for each ^4 e A, S~^AS is a diagonal 
matrix. Therefore A is semisimple and has a second basis Eo,...,Ed (see [23]). These are 
matrices satisfying 

Eo^-J 
n 

EiEj — SijEi 



T.E^ = I. (2-7) 

j=0 

The matrix ^ J (where J is the all-one matrix in A) is a minimal idempotent (idempotent is 
clear, and minimal follows from the rank (J =1)). The Ei, for (0 < i, j < d) are known as the 
primitive idempotent of Y . Let P and Q be the matrices relating our two bases for A: 

d 

Aj = J2PijEi, 0<j<d, 

i=0 



continuous-time Quantum walk 10 



1 

Ej = -Y.QijA, 0<j<d. (2-8) 



Then clearly 



It also follows that 



PQ^QP^ nl. (2-9) 



AjEi = PijEi, (2-10) 

which shows that the Pij (resp. Qij) is the i-th eigenvalue(resp. the i-th dual eigenvalue ) 
of Aj (resp. Ej) and that the columns of Ei are the corresponding eigenvectors. Thus = 
rank(£^j) is the multiplicity of the eigenvalue Pij of Aj (provided that ^ Pkj for k ^ i). 
We see that mo = 1, Y^i^^i — and mj =trace£'j = n{Ei)jj (indeed, E^ has only eigenvalues 

and 1, so rank(_E'fc) equals the sum of the eigenvalues). Also, by [25, 26], the eigenvalues and 
dual eigenvalues satisfy 

PiO — QiO — 1) Poi — ki, Qoi — TUi 

rrijPji = kiQij, <i,j <d, (2-11) 
where for all integer i {0 < i < d), set = and note that k^ ^ 0, since Ri is non-empty. 

2.3 The Terwilliger algebra and its modules 

We now recall the dual Bose-Mesner algebra of Y. Given a base vertex a eV, for all integers 

1 define E* — E*{a) e Matv{C) {0 < i < d) to be the diagonal matrix with {(3, /3)-entry 

r 1 if {a,p) e Ri, 
= (« e V). (2-12) 

i otherwise 

The matrix E* is called the i-th dual idempotent of Y with respect to a. We shall always set 
= for i < or i > (i. From the definition, the dual idempotents satisfy the relations 

d 
i=0 
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E*E] = SijE* 0<t,j<d. (2-13) 

It follows that the matrices Eq, E^, .... form a basis for the subalgebra A* = A* (a) of 

Matv{R)- A* is known as the dual Bose-Mesner algebra of Y with respect to a. For each 

integer i {0 < i < d),let A* — A* (a) denote the diagonal matrix in Matv{R) with /3)-entry 

{A^)^p,^)^n{E,)^,^ {(3eV). (2-14) 

With reference to [25, 26] the matrices A^, A^, A^ form a second basis for A* and satisfy 

A*o = I, 
Af = At, 
Al + A* + ...+A*^nE*, 

AiA- = E (2-15) 

h=0 

Then by combining (2-8) with (2-12) and (2-14) we have 

d 

A] = Y,Q,^El 0<j<d, 



i=0 

d 



E*^-Y,PioAl 0<J<d. (2-16) 

^ i=0 

Let Y = (y, {-Ri}o<i<d) denote a scheme. Fix any a ^ V , and write A* = A*(«). Let T = T{a) 
denote the subalgebra of Matv{C) generated by A and A*. We call T the Terwilliger algebra 
of Y with respect to a. 

Thus, we can define quantum decomposition for distance regular graphs by the following 
lemma: 

Lemma (Terwilliger [17]). Let F denote a distance regular graph with diameter d. Fix any 
vertex a of F, and write E* = E*{a) {0 < i < d), Ai = A andT = T{a). Define A" = A- (a), 
= AO(q;), A+^A+{a) by 

A- = Y: EUAEl A' = Y. EtAEl = ^ E*^,AE*. (2-17) 

i=l i=l i=l 
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Then 

A^A+ + A- + A^, (2-18) 
where, this is the quantum decomposition of adjacency matrix A such that, 

{A-y - A+, {Ay = (2-19) 

which can be verified easily. 

Let W — denote the vector space over C consisting of column vectors whose coordinates 
are indexed by V and whose entries are in C. We observe Maty{C) which acts on W by left 
multiplication. We endow W with the Hermitian inner product (, ) which satisfies (u, v) — u^v 
for all rt, V e , where t denotes the transpose and - denotes the complex conjugation. For 
all P & V, let 1/3) denote the element of W with a 1 in the (3 coordinate and in all other 
coordinates. We observe {|/3)|/3 £ V} is an orthonormal basis for W. Using (2-12) we have 

E*W = span{\(3)\(3 e V, (a, (3) e Ri}, < i < d. (2-20) 

With the use of Eq.(2-20) and since £ V} is an orthonormal basis for W, we get 

w = E*w e E^w e • • • e e^w, (2-21) 

where the orthogonal direct sum, with < i < d, E* acts on W as the projection onto E*W, 
similarly W = Yh=o EiW . We call E^W the z-th sub constituent of T = (V, R) (distance 
regular graph) with respect to a. For < i < o? we define 

\<I>^)-Y.\(^). (2-22) 

where the sum is over all vertices P such that (a, P) e Ri. We observe e EqW. 

By a T-module we mean a subspace U CW such that TU C U. Let U denote a T-module. 
Then U is said to be irreducible whenever U is nonzero and U contains no T-modules other 
than and U. Let U denote an irreducible T-module. Then U is the orthogonal direct sum 
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of the nonzero spaces among E^U, E^U, E^U ([17], Lemma 3.4). By the endpoint of U we 
mean min {i|0 <i <d, E*U ^ 0}. By the diameter of U we mean |{i|0 < i < E*U 7^ 0}| - 1. 
We say V is thin whenever E\U has dimension at most 1 for < i < d. There exists a unique 
irreducible T- module which has endpoint ([34], Prop. 8.4). This module is called Wq. For 
< i < d the vector 10^) of Eq.(2-22) is a basis for E\Wq ([17], Lemma 3.6). Therefore Wq 
is thin with diameter d such that the module is orthogonal to each irreducible T-module 
other than Wq ([35], Lem. 3.3). 

2.4 Stratification 

For q; 7^ /3 let 9(q;, /3) be the length of the shortest walk connecting a and /3. By definition 
9(q;, /3) = for all a E V. The graph becomes a metric space with the distance function d. 
Note that d{a,P) — 1 if and only if a ~ /3. We fix a point o e F as an origin of the graph. 
Then, the graph is stratified into a disjoint union of associate classes rfc(o): 

y = U Fkio), Fkio) ^{aeV; d{o, a) = k}. (2-23) 

fe=0 

With each associate class rfe(o) we associate a unit vector in P{V) defined by 



where, \a) denotes the eigenket of a-th vertex at the associate class Tkio) and = |rfc(o)|. 
Obviously, the two sites, first and last strata, have the same stratification. The closed subspace 
of l'^{V) spanned by is denoted by A(G). Since {\4>k)} becomes a complete orthonormal 

basis of A(G), we often write 



Let Ai be the adjacency matrix of a distance regular graph F = {V,R) for vacuum state |0) 




(2-24) 



A(G) = E©C|0,). 



(2-25) 



k 



we have 



Ak\0)= E 1/5)- 



(2-26) 



/36rfc(o) 
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Also, we have 



A,\0) e EtW 



(2-27) 



Then by using unit vectors (|0) = |o) = \(po), with o E V as the fixed origin), and 
Eq.(2-26), (2-27) we have 



3 Continuous-time quantum walk on graphs with known 
spectrum 

The continuous-time quantum walk is defined by replacing Kolmogorov's equation with Schro- 
dinger's equation. Continuous-time quantum walk was introduced by Farhi and Gutmann [5] 
(see also [8, 30]). Our treatment, though, follow closely the analysis of Moore and Russell 
[30] which we review next. Let 1^{V) denote the Hilbcrt space of C-valued square-summable 
functions on V. With each a &V we associate a ket defined by \a), then {|a), a e V} becomes 
a complete orthonormal basis oil'^{V). 

For <i < d the vector of Eq.(2-24) is a basis of E*Wo, where Wq is unique irreducible 
T- module which has endpoint ([34], Prop. 8.4). Therefore, Hilbert space of continuous- 
time quantum walk starting from a given site corresponds to the irreducible ( walk starting 
site related T-algebra) T- module Wq with maximal dimension. Hence other irreducible T- 
modules of Terwilliger algebra T are orthogonal to Hilbert space of the walk. 

Let |0(t)) be a time-dependent amplitude of the quantum process on graph F. The wave 
evolution of the quantum walk is 



Ak\(f>o) 



(2-28) 



ih-m)) = H\m, 



(3-29) 
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where we assume h = 1, and \4>q) is the initial amphtude wave function of the particle. The 
solution is given by |0o(O) — ^~^^*\4>o)- It is niore natural to deal with the Laplacian of the 
graph, defined as L — A — D, where D is a diagonal matrix with entries Djj = deg{aj). This 
is because we can view L as the generator matrix that describes an exponential distribution of 
waiting times at each vertex. But on s-regular graphs, D — ^I, and since A and D commute, 
we get 

This introduces an irrelevant phase factor in the wave evolution. In this paper we consider 
L — A — Ai. Then using Eq.(2-8) we have 

\Mt)) = e-'^'\<Po) = e-^i:to^"^>*|0o), (3-31) 

where using the algebra of idempotents i.e., Eq.(2-7), the above amplitude of wave function 
can be written as 

l0o(t))=Ee-^^-*£;#o). (3-32) 

Now using Eqs.(2-8), (2-9), (2-24) and (2-28), the matix elements of idempotent operators 
between eigenstates strata and eigenstates of vertices can be calculated as 

{<t>k\Ei I 0o) = {<l>k\-^QiiAi I 0o) = - ^QiiiM^i I 0o) = ^^Qki, (3-33) 



n n n 



or 



im I <Po) = -Qki, (3-34) 
n 

for every e Tk{o). 

Finally multiplying (3-35) by and using (3-38) we get the following expression for the 
amplitude of observing the particle at vertex (3 at time t 

iPlMt)) = Ee-^-*(/3|£;, I 0„) = i^e-^"*g,,. (3-35) 

1=0 i=0 

Obviously the above result indicates that the amplitudes of observing walk on vertices belong- 
ing to a given stratum are the same. Actually one can straightforwardly deduce from formula 
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(3-38) that the transition probabihties between the vertices depend only on the distance be- 
tween the vertices (kind of relations or association classes) , irrespective of which site the walk 
has started. This is due to the fact that in association schemes the coloring of underlying 
graphs or the set of relations between vertices determine thoroughly everything. Hence for a 
continuous-time quantum walk over a graph associated with a given scheme with diameter d, 
we have exactly [d+ 1) different transition probabilities (i.e., the number of strata or number 
of distinct eigenvalues of adjacency matrix). 

At the end by straightforward calculation, one can evaluate the average probability for 
finite graphs of association schemes as 



for every \p) e Tk{o). 

4 Continuous-time quantum walk on group schemes 

In this section we briefly discuss continuous-time quantum walk on group schemes with real 
and complex representations separately. 

4.1 Group association schemes 

In order to study the continuous-time quantum walk on group graphs, we need to study the 
group association schemes. One of the most important sources of association schemes are 
groups. Let G be a group acting on a finite set V. Then G has a natural action on V x V 
given by g{a, (3) = {ga, gP) ior g e G and a,P eV. The orbits 



oiVxV are called orbitals. The group G is said to act generously transitive when for every 
pair {a, (3) e V xV there is a group element g & G that exchanges a and (3, that is ga — (3 and 




(3-36) 



{{ga,gP)\geG} 



(4-37) 
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gP = a. When G acts generously transitive, the orbitals form the relations of that association 
scheme. Now, in the following, we consider the orbitals which correspond to the conjugacy 
classes of G. Let G be a finite group, Cq = {e}, Ci, Cd the conjugacy classes of G. Let GxG 
act on G with the action defined by P{ai, a2) — a^^/9Q!2 where ai, 0:2 G G. Then the diagonal 
action of GxG on GxG is given by 7)(q;i, — (/3(q;i, q;2), 7(0:1, a2)) — {ai^Pa2, ^7q;2). 
One can show that {j3i, (32), (71, 72) G G x G belong to the same orbital of G x G if and only if 
Pi^P2, li^l2 belong to the same conjugacy class of G. Thus in this case the orbitals correspond 
to the conjugacy classes of G. For i = 0, 1, d define 

Ri^{{a,P)\a-'peCi}, (4-38) 

then Ri are the orbitals of G x G and hence X(G) — (G, {Ri}o<i<d) becomes a commutative 
association scheme and it is called the group association scheme of the finite group G [26] . We 
define class sum Ci ior i — 0,1, ...,d as 

G, = ^ 7 e GG, (4-39) 

then, for regular representation we have Ci\a) — J^-yeCi Therefore in regular representa- 
tion, the classes sum Ci{i = 0,1, ...,d) have the following matrix elements 

r 1 if (a,/3) eG,, 
(a)a,/3 = {a,PeG). (4-40) 

I otherwise 

Comparing the above matrix elements with those of adjacency matrices given in (2-4), we see 
that the classes sum are the corresponding adjacency matrices of group association scheme 
with the relation defined through conjugation. It is well known that the classes sum of finite 
group G form the basis of center of its GG ring which is certainly a commutative algebra, 
hence they are closed under multiplication defined in GG, i.e, we have 

k=0 
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(see details in [28]), where p^j {i,j,k = 0,1, ...,0?) are intersection numbers of the group 
association scheme X{G) and have the following form: 



k _ \Ci\\Cj\ ^ x{ai)x{aj)xM . 
\G\ ' ^ ^ 

where the sum is over all the irreducible characters x G [29]. Therefore, the idempotents 
{Eq, El, Ed} of the group association scheme X{G) are the projection operators of GG- 
module i.e, 

Ek = ^j: Xk{a-')a. (4-42) 
Thus eigenvalues of adjacency matrices of Ak, and idempotents Ek, respectively are 

Pik = ^^XiM: Qik = dkXk{0ii): (4-43) 

rrii 

where dj = Xj(l). The above defined group scheme is in general non-symmetric scheme and 
it can be symmetric provided that we choose a group whose whole irreducible representations 
of chosen group are real, such as symmetric group S^. In appendix A, we have explained how 
to construct a symmetric group association scheme from a non-symmetric one. 

A. Continuous-time walk on group schemes with real representations 
In a finite group G with real conjugacy classes Co — {e},Ci, ...,Crf, i.e., G{a) — G{a~^) for 
all q; e G, all irreducible characters Xi ^-^e real. Thus using (4-43) we can study continuous- 
time quantum walk on its underlying graph, where the amplitude of observing the particle at 
stratum k at time t, i.e, Eq.(3-35) reduces to 

{<t>k\Mt)) = ^Ed.e""'^"/'"'*3^, (4-44) 
also, the average probabilities over large times becomes 

P{k) = ^Y.d'Mak)\\ k = 0,l,...,d. (4-45) 

^ i=Q 

Therefore, the probability of observing the walk at starting vertex, i.e., the staying probability 
is 

m = 5E^?lx.(0)r = ;^E4 (4-46) 
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As examples we will study continuous-time quantum on G — Sn, graphs in section 6. 
B. Continuous-time walk on group schemes with complex representations 

In general all conjugacy classes of a given finite group are not real, hence some of its irreducible 
representations become complex and consequently we encounter with directed undelying graph 
or non-symmetric association scheme. But following instruction of appendix A we can generate 
a symmetric association scheme out of non-symmetric association scheme. Thus in this case, 
for continuous-time quantum walk, we need to use formulas (A-v) and (A-vii) of the appendix 
A, where the amplitude of observing the particle at stratum k at time i.e., Eq.(3-35) is 



/at ; -idj)cixi(ai)t ^ 

■^Si=o^i6 XiyO'-k) for real representation 

^y— d+l ^ -id,;)ti(x,:(ai)+X,:(ni))t 



^ die {Xi{cik) + Xi(afc)) for non-real representation. 

(4-47) 



Also, the average probabilities are 



~ I]i=o d}\Xi{oik)\'^ for real representation 



d+l 



^ 0?- |(Xi(Qife) + Xi{oik))?- for complex representation. 

In this case the staying probability is 



(4-48) 



d+l d+l 

„ I 2 1(2 



^(0) = ;|(E^.lx.(0)r + 4 E dlUrm") = -,{Y.dt + ^ E <)■ (4-49) 

" i=0 1=1+1 " 1=0 i=l+l 

As an example we will study G = C„ in section 6. 



5 Investigation of Continuous-time quantum walk on 
distance regular graphs via spectral distribution of 
adjacency matrix 

Even though in general modules of TerwiUiger algebra stratifies the underling graph of cor- 
responding association scheme but in case of distance regular graphs, due to the existence of 
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raising and lowering operators given by (2-17), we can have a stratification similar to the one 
introduced in Ref[l]. As we will see in the following this particular kind of stratification leads 
to introduction of spectral distribution for adjacency matrix. 

Here in this case, the strata states \(f)k) are the same as defined by (2-24) of subsection 2.4, 
but further using Eqs.(2-6) and (2-18) one can show that, the raising and lowering operators 
given by (2-17) act over them as follows 

1 = <Pk+i), k>0 (5-50) 

A-\cf>o)^0, A-\cj>k)^V^\cj>k-i), k>l (5-51) 
AVfe) = («fc+i)|0fe), k>0. (5-52) 

As mentioned in section 2, \(f)k), k — 0,1, form basis for Wq which is the irreducible 
T-module with maximal dimension, therefore all basis of the irreducible T-module Wq can 
be obtained by repeated action of raising operator A'^ on reference state |0o) ) and we have 
cuk — Ck-ibk, CKfc+i = ai + Ck + bk similar to Ref.[l]. The space Wq equipped with set of 
operators (T , A'^ , A~ , A'^) is an interacting Fock space associated with the Jacobi sequence 
{ck-ibk,ai-bk- Ck, A; = l,2, ...}. 

Usually the spectral properties of the adjacency matrix of a graph play an important role 
in many branches of mathematics and physics and the spectral distribution can be generalized 
in various ways. In this work, following Ref.[22], the spectral distribution fj, of the adjacency 
matrix A (where A = Ai) is defined as: 

< A"^ >^ I x'^nidx), m = 0,l,2,... (5-53) 
Jr 

where < . > is the mean value with respect to a the ground state |0o) , and according to Ref. [22] , 
the < A^ > coincides with the number of m-step walks starting and terminating at a. Then the 
existence of a spectral distribution satisfying (5-53) is a consequence of Hamburgers theorem, 
see e.g., Shohat and Tamarkin [[38], Theorem 1.2]. 
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We may apply the canonical isomorphism from the interacting Fock space (Hilbert space of 
continuous-time quantum walk starting from a given site, i.e., strata states or more precisely 
Wq the irreducible T- module with maximal dimension starting) onto the closed linear span 
of the orthogonal polynomials determined by the Szego-Jacobi sequences {{u!k}, {(^k}), where 
for distance-regular graphs, the parameters cuk and ctfe are defined as 

<^ik = Cjfc_i6jfc, Qijfc = ai - bk-i - Cjk_i, (5-54) 

where ak,bk and Ck have been introduced in the relations (2-3). More precisely, the spectral 
distribution /i under question is characterized by orthogonalizing the polynomials {Qn} defined 
recursively by 

xQnix) = Qn+lix) + an+lQnix) + UnQn-li.x), (5-55) 

for n > 1. If such a spectral distribution is unique (for distance- regular graphs is possible), 
the spectral distribution fj, is determined by the identity: 

= L^^ .^fiW.t^. (5-56) 

jRX-y a: - Qi - ^ ' Qn{x) ^^x-Xi 



^3 x — ceA — 



where Gn{x) is called the Stieltjes transform and Bi is the coefficient in the Gauss quadrature 
formula corresponding to the roots xi of the polynomial Qn{x) and where the polynomials 
{Q^n^} are defined recursively as 

Q^\x) = 1, 

Q^i\x) ^x - a2, 

xQm{x) = Q^^Ux) + a„+2Q«(^) + a;„+itf2i(x), 
for n > 1. 

Now if G^{x) is known, then the spectral distribution can be recovered from G^{x) by 
means of the Stieltjes inversion formula: 

1 fv 

fj,{y) — fj,{x) — lim / Im{GJu-\-iv)}du. (5-57) 

TT V — >0+ Jx 
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Substituting the right hand side of (5-56) in (5-57), the spectral distribution can be determined 
in terms oi xi,l — 1,2, the roots of the polynomial Qn{x), and Guass quadrature constant 
Bi, 1 — 1,2, ... as 

li^Y.BiS{x-xi) (5-58) 

(for more details see Refs.[36, 38, 39, 40].) 

Finally, using the relations(2-28) and orthogonal polynomial Pn{x), where they satisfy the 
recursion relations (2-6) and (5-55), the other matrix elements {(l)k\A^ \ 0o) can be written as 

I 0o) = ^ /■ x"'Pk{x)pL{dx), m = 0, 1, 2, .... (5-59) 



/ak 

Our main goal in this paper is the evaluation of amplitude for continuous-time quantum walk 
by using Eq.(5-59) such that we have 

{(PklMt)) = ^ / e-'^'Pk{x)ii{dx), (5-60) 



where (0fc|0o(^)) is the amplitude of observing the particle at level k at time t. The conservation 
of probabihty X^fe=o I (0fe|0o(^)) 1^ = 1 follows immediately from Eq.(5-60) by using the com- 
pleteness relation of orthogonal polynomials P„(a;). Obviously evaluation of (0fc|0o(^)) leads 
to the determination of the amplitudes at sites belonging to the associate scheme (stratum) 
Tk{o). As proved in the appendix B, the walk has the same amplitude at all sites belonging 
to the same associated class ( stratum). Also, for the finite graphs, the formula (5-60) yields 

(0fc|0o(t)) ^^Y.Bie-''^''Pk{xi), (5-61) 

where by straightforward calculation one can evaluate the average probabilities for the finite 
graphs which yields 

m = lim ^ r I {<Pk\Mt)) r dt = - E^fn'(^0- (5-62) 

T-»oo I Jo ttk 1 
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6 Examples 



Here in subsection 6.1 we study continuous-time quantum walk on graphs with known spec- 
trum. 



6.1 Complete graph Kn 

A complete graph with n vertices (denoted by Kn) is a graph with n vertices in which each 
vertex is connected to the others (with one edge between each pair of vertices). A complete 
graph is the trivial association scheme, with the intersection numbers 



ai — n — 1, bi — 1, Co — n — 1. 



(6-63) 



It is straightforward to show that its adjacency matrix A has egivenvalues 1 and —1 with 

degeneracies 

Eo — —J, 

n 



n 



n-1 -1 
-1 n-1 



-1 
-1 



-1 ■■■n-1 



(6-64) 



and using Eq.(3-35) we have 



^(e-** + {n- l)e^) for /3 = 



-2 



-it(n-2) 



(6-65) 



■ sm{^)c^(^ for/3^0. 
Finally using Eq.(3-36) we get the following expressions for the average probabilities 



1 - 



2(n-l) 



for /? = 
for /3 ^ 0. 



(6-66) 
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6.1.1 Symmetric group Sn 

The symmetric group Sn is ambivalent in the sense that C{a) = C{a~^) for all a e Sn, 
therefore its conjugacy classes form a symmetric association scheme (see Fig.2). 

For group Sn, conjugacy classes are determined by the cycle structures of elements when 
they are expressed in the usual cycle notation. The useful notation for describing the cycle 
structure is the cycle type [i/i, i/2, i^n] , which is the listing of number of cycles of each length 
(i.e, vi is the number of one cycles, 1^2 is that of two cycles and so on). Thus, the number of 
elements in a conjugacy class or stratum is given by 

'^1-— ]' = .,!2-.2!.'..n-.„!- ^^-^^^ 
On the other hand a partition A of n is a sequence (Ai, A„) where Ai > • • • > A„ and 
Ai + • • • + A„ = n, where in terms of cycle types 

Xl = Ui + U2-\ ^Un, 

Al = 1^2 + i^S H \-l^n, 



Xn = l^n- (6-68) 

The notation X\- n indicates that A is a partition of n. There is one conjugacy class for each 
partition A h n in 5"^, which consists of those permutations having cycle structure described 
by A. We denote by Cx the conjugacy class of Sn consisting of all permutations having cycle 
structure described by A. Therefore the number of conjugacy classes of Sn, namely diameter of 
its scheme is equal to the number of partitions of n, which grows approximately by 4^^ ^""^^"^^- 
We consider the case where the generating set consists of the set of all transposition, i.e, 
Ci = C[2,i,i,i,i...,i]- For the characters at the transposition, it is known that [31] 



2\{n - 2)\dim{px) ^ 













V 


y 2 ) 




[ 2 J 





(6-69) 
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Here, A' is the partition generated by transposing the Young diagram of A, while A^- and Xj 
are the j-th components of the partitions A' and A, and px is the irreducible representation 
corresponding to partition A. 

Then the eigenvalues of its adjacency matrix can be written as 



Al 



dxh, 
mx ' 





( 








\ 


E 












3 


\ 






[ 2 J 


/ 



(6-70) 



Therefore by using Eqs.(4-44) and (6-70), one can obtain the amplitudes on symmetric groups. 
As an example we obtain amplitude for associate class of conjugacy class of n-cycles, as 

(2i sin(nt/2))"-i 



{(t>n\Mt)) = 



'nn\ 



(6-71) 



where the results thus obtained are in agreement with those of Ref [32] . 

In the above calculation, we have used the following results for the characters of the n-cycles 



Xx{{n)) 



_lY-k for A = (A;, 1, . . . , 1), A; e {1, . . . , n} 
otherwise 



and 



n — 1 



k-l 



Pxi = -{2nk — — n). 
2 



X(fc,i,-,i) (id) = rf^m(p(fc,i,.,i)) = 
6.1.2 Dihedral group 

The dihedral group G — is semi-direct product of cyclic groups and Z2 with cor- 
responding generators a and h. Hence it is generated by generators a and h with following 
relations: 

D2m = (a, 6 : a"' = 6^ = 1, h'^ah = a"^). (6-72) 

In finding its conjugacy classes, it is convenient to consider whether m is odd or even. 
Hence, we will study continuous-time quantum walk on dihedral group for odd and even m, 
separately. 
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1. m=odd the dihedral group has precisely ^(m + 1) conjugacy classes: 



Co = {1}, Ci = {6, ab, a% a'^-^h], 
C2 = {a, a-^}, = {a("^-^)/^ a-("^-^)/2| 

for more details see Ref.[28]). 
By using Eqs.(4-43) and (3-35) we obtain the following amplitudes 



^((m - 1) + cos(mt)) 
^(-isin(mi)) 



(6-73) 



for A; = 

for A; = 1 (6-74) 

^(cos(mt) - 1) for A; = 2, 3, (m + 3)/2. 

Uusing Eq.(3-36) we obtain the following average transition probabilities to strata k at time t 

(6-75) 



P{k) = <^ 



1 

2m 



for A; = 
for A; = 1 
for A; = 2,3, ...,(m + 3)/2, 



1. m=even the dihedral group {iti — 21) has precisely I + 3 conjugacy classes: 
Co = {1}, Ci = {a'}, C2 = {a, a''}, Q = {a'-^'«"'^'} 



Ci+i = {a'^ft : < J < Z - 1}, Ci+2 = {a'^+'& : < j < Z - 1} 



(6-76) 



( for more details see Ref.[28]). Since we will restrict our attention to Ci generate group, 
therefore we consider conjugacy classes as 

Co — Co, Ci = C;+i U C;+2, C2 = Ci, C3 = C2, C4 = C3, • • • , C;+i = Q. (6-77) 

In this case, calculation of amplitudes and probabilities is similar to that of dihedral groups 
with odd m. 
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6.1.3 Cycle graph C„ 

A cycle graph or cycle is a graph that consists of some number of vertices connected in a closed 
chain. The cycle graph with n vertices is denoted by C„, where its graphical representation 
cyclic group Zn — (a), whit a"' — 1. In this case, we consider the orbitals to correspond 
to the conjugacy classes of cyclic group. Also, we give C = C{a) U C{a~^) for all a G Zn, 
therefore the relations Ri form a symmetric association scheme with d classes on C„ (called 
the conjugacy scheme of C„). Let cvj — e^'^'^jh for j = 0, 1, n — 1. Using the properties of 
characters of cyclic group, and Eq.(A-iv) and (A-v), we obtain 

^ii = Xi(l) + Xi(^ - 1) = ^3 + ^T^ = 2 cos(27rj7n). 

Qkj = XkU) + Xk{n -j) = 2 cos{2ti jk/n). (6-78) 
Therefore by using Eq.(3-35), for strata k we have 

^ (e-^* + 2 Ef-n e-^*'^°'(2''^/")) for A; = 

{MMt)) = \ " ' (6-79) 

[ ^(e"'* + 2 E?=o e-^*^°«(2^^'/") cos(27rjfc/n)) for = 1, 2, d, 

where the results thus obtained are in agreement with those of Ref.[l, 33]. Thus, one can 
evaluate the average probability of staying at origin for large time as follows Eq.(4-49) 

P(0) = l(l + 4X:cos^(0))) = ^(l + 4d), (6-80) 

and using Eq.(3-36) we get the following probabilities of transition to stratum k 

P{k) = ^(l + 4Vcos'(27rjA;/^))). (6-81) 
^ j=i 

Where here we have considered odd, n, and calculation for even n is similar to that of cycle 
graph with the odd one. 

In the remaining the part of paper we will consider the examples of continuous-time quan- 
tum walk which can be investigated via spectral distributions of their corresponding adjacency 
matrices. 
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6.2 Strongly regular graphs 

A graph (simple, undirected and loopless) of order n is strongly regular with parameters 
n, A, ?7 whenever it is not complete or edgeless and 

(i) each vertex is adjacent to k vertices, 

(ii) for each pair of adjacent vertices there are A vertices adjacent to both, 

(iii) for each pair of non-adjacent vertices there are 77 vertices adjacent to both. 
For strongly regular graph, the intersection numbers are given by 

ai — k; bi — 1, b2 — rj; cq — k, ci — k — X — 1. (6-82) 

By using formula (5-58) one can straightforwardly get the following spectral distribution 

/i = Bi6{x - xi) + B26{x - X2) + Bs6{x - x^), (6-83) 

where, we obtain Xi and for i = 1, 2, 3 respectively as 

Xi — K, 



X3 = ^(A - 77 - V(A- 77)2- 4(77 -«)), (6-84) 



— k{\ — rj) + {rj — kY 



^ -K^^{X - rjf - 4{rj - k) + k{X - rj) + 2k 

^ (A - 77 - 2k)^{X - rj)^ - 4(77 -k) + {X- rj)^ - 4(77 - k) ' 

i^\/(>< - TlY - - «) + -r]) + 2K 

S3 = ^ , . (6-85) 

(-A + ?7 + 2k)J{X - rj)^ - 4(7/ - /t) + (A - 7/)2 - 4(7/ - k) 



Again using Eq.(5-60) one can obtain the amplitudes for quantum walk at strata k and time 
t. 

For example we study continuous-time quantum walk on the following graphs. 
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6.2.1 Petersen graph 

Petersen graph is strongly regular graph with parameters (n, k, A, 77) = (10, 3, 0, 1) (see Fig.3). 
The intersection numbers and spectral distribution are 

oi =3, 02 = 6; 61 = 62 = 1; Co = 3, ci = 2. 

= T^^i^ - 3) + - 1) + + 2). (6-86) 
iU z o 

Therefore, the amplitudes for walk at time t are 

(02|0o(t)) = ^ \y - 3)e-'^{dx) = -^i-e-'' + ^e^^* + ^e^^*). (6-87) 

6.2.2 Normal subgroup graphs 

Let G be a finite group, and P = {Pq, Pi, Pd} be a blueprint of it. We always assume that 
the sets Pj are so numbered that the identity element e of G belongs to Pq. If Pq = {e}, 
then P is called homogeneous. Let {Rq, Ri, Rd} be the set of relations Ri = {(a,/?) e 
G<^G\a~^P e Pi} on G. Now, we define a blueprint for group G which form a strongly regular 
graph. If if is a subgroup of G, we define the blueprints by 

Po = {e}, Pi^G-{H}, P2^H-{e}. (6-88) 

This blueprint form a strongly regular graph with parameters (n, k, A, 77) = (|G|, |G| — |G| — 
2\H\,\G\ - \H\). 

As an example, we consider G — D2m'- 

1. m=odd in this case the subgroup is defined as 

H = {e, a, a-\ a^^-^^'\ a'^^''-^/''}. (6-89) 
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Therefore the blueprints are given by 

Po = {e}, Pi = {b,ab,a\...,a"'-^b}, 
P2 = {a, a-\ a(™-^)/^ a-^'""^)/"}, (6-90) 

which form a strongly regular graph with parameters (2m, m, 0, m). By using Eqs. (6-82), (6-83), 
(6-84) and (6-85), we get the following expressions for the intersection numbers and spectral 
distribution 



Qi — m, a2 — m — 1; bi — 1, b^ — m; Cq — m, Ci — m — 1. 

If/ N ~ 1 r/ \ 1 c-/ \ 

A* = 7: — 0{x - mj H d(x) + - — d{x + m). 

2m m 2m 

Therefore, the amplitudes for walk at time t are 

(0o|0o(^)) = — ((m - 1) + cos(mi)) 



m 



{4>i\Mt)) 



'—i sin(mt)) 



m 



cos(mt) — 1). 



Also using Eq.(5-62) we evaluate the average probabihties as 

^((m-l)2 + i) forA; = e 

m = 



1 

2m 

3(m-l) 
2m2 



for /c = 1 
for k^2. 



2. m=even we consider m = 21 and define subgroup as 



H = {e, a, a 



a'}. 



Therefore the blueprints are given by 



Po = {e}, Pi = {a2^6,a2J+^6, < j < Z - 1}, 



P2 = {a, a-\ a('-^), a-^'"^), a'}. 



(6-91) 



(6-92) 



(6-93) 



(6-94) 



(6-95) 



which are strongly regular graphs with parameters (2m, m, 0, m). In this case, the calculation 
is similar to that of dihedral groups with m, odd. 
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6.3 Cycle graph Cn 

spectral distribution 

A cycle graph or cycle has already been defined in the subsection 6.1. The cycle graph C„ 
is the distance regular graph, and quantum walk on them turns out to be different for odd and 
even n, hence, below we treat them separately. 

Odd n. For odd n — 2d+l, the diameter of graph C2d+i is d also the intersection numbers 

as 

oi = 02 = • • • = = 2; 6i = 62 = • • • = &d = 1; Co = 2, ci = C2 = • • • = Cd-i = 1. (6-96) 
In this case, the adjacency matrices are given by 

Ao = In: Ai = S^ + S-\ A2 = S^ + S-\---,Ad = S'' + S-^ (6-97) 

where S is the shift operator 5'2<^+i = 1. Now by using the Bose-Mesner algebra of distance- 
regular graph for adjacency matrices (2-6), we obtain 

Ai = 2Ti{A,/2), />1, (6-98) 

where 7] is Tchebichef polynomials of first kind. Also, Eq.(6-98) and Eq.(2-6) ior i — d we 
have 

Td{Ai/2)^Ta+i{A,/2). (6-99) 

The canonical isomorphism from the interacting Fock space onto the closed linear span of the 
orthogonal polynomials maps states \(f)i) into orthogonal polynomials Pi{x) — 2Ti{x/2). 
It is straightforward to show that the spectral distribution Eq.(5-58) as 

li=— ^5(x- 2) + ^^ yS(x- 2 cos(^^)) (6-100) 

where xi — 2cos(^^) are the roots of Eq.(6-99). Therefore, the amplitudes for observing the 



walk at time t and /cth associated class is 

(0fc|0o(t)) = ^ / e--*/2r,(x/2)/i(dx) 

\ dh JR 
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= -^{e-'' + 2j2e-''"''^^''^^^'^^'^ cos{2kln/i2d + 1))), (6-101) 

2(2 + 1 

where the results thus obtained are in agreement with those of Refs.[l, 33]. 

Even n. For even n — 2d, the diameter of graph C2d is d and the intersection numbers are 

oi = 02 = • • • = a<i_i = 2, = 1, 6i = 62 = • • • = bd-i = 1, 6^ = 2, 

Co = 2, ci = C2 = • • • = Cd_i = 1. (6-102) 

In this case, its calculation is similar to that of cycle graph with n, odd. 

In the limit of large n, the cycle graph C„ is the same as the infinite line graph, with the 
intersection numbers 

ai = 02 = • • • = 2; fei = 62 = • • • = 1; Co = 2, ci = C2 = • • • = 1. (6-103) 
Then the orthogonal polynomials are Pi{x) — 2Ti{x/2) and the spectral distribution n is 

pt^-—L=, -2<x<2,. (6-104) 

TT V4 — X'' 

Then one can obtain spectral distribution and amplitude for every vertex at the time t, where 
the results thus obtained are in agreement with those of Ref.[l]. 

6.4 Johnson graph 

Let v,dhe a pair of positive integers such that d < v. Put -S" = {1,2, ...,v} and V = {x G 
S :| X \— d}. We say that x,y & V are adjacent if d— \ x Hy \= 1. Thus a graph structure 
is introduced in V, which is called a Johnson graph and denoted by J{v,d) (see Fig.4). By 
symmetry we may assume that 2d < v. Consider the growing family of Johnson graphs J(f , d), 
where 0? — > 00 and ^ ^ p e (0, 1] . Then the associated orthogonal polynomials are (for more 
details see Ref.[36]) 

A. for p = 1, we have Laguerre polynomials Ln{x) with the following recurrence formula: 

Lq{x) = 1, 
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Li{x) = X — 1, 

) + (2n + l)Ln{x) + n'^Ln_i{x), n > 1. (6-105) 

By using the fact that the Laguerre polynomials are orthogonal polynomials with respect to 
the spectral distribution e~^dx and following paper [1] we obtain the following amplitudes 



i<f>Mt)) - (6-106) 

B. for < p < 1, by modifying the Meixner polynomials Mn{x), we have the recurrence 
formula: 

Mo{x) - 1, 

Mi{x) = X, 
2n 

xMn{x) = Mn+i{x) + . Mnjx) + n^M„_i(x), n>l. (6-107) 

- p) 

Hence by using the fact that the Meixner polynomials are orthogonal polynomials with respect 
to the spectral distribution E^o '^^^{^YK^ - =£±B=^) and Eq.(5-60) one can obtain 

P P p(2 p) 

the amplitudes. As an example, we obtain the amplitude at the origin at time t as 

2(1 -p), p 



(0o|0o(t)) = E^3f^(^)'e" ^ ■ (6-108) 

f^^Q Z p Z p 

6.5 Product of association schemes 

In this section, we recall some basic facts about the symmetric product of trivial schemes, (see 
[37] for more details.) This product is important not only as a means of constructing new 
association schemes from the old ones, but also for describing the structure of certain schemes 
in term of particular sub-schemes or schemes whose structure may already be known. Then 
using Eq.(3-35) and (5-60), we can evaluate amplitudes of continuous-time quantum walk 
on new association schemes. The symmetric product of o?-tuples of trivial scheme Kn with 
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adjacency matrices of In, Jn — In is association scheme with the following adjacency matrices 
( generators of its Bose-Mesner algebra) 



permutation 



Ai^ ^ (Jn - In) ® {Jn - Q-- ® {Jn - In) ®In ® ■■■ ® 4, (6-109) 
permutation ". 

where J„ is n x n matrix with all matrix elements equal to one. This scheme is the well known 
Hamming scheme with intersection number 

^_^(n-lMd-l)...(^-. + l)_ ^^^^^ 
i! 

bi — i, 1 <i < d, 

Q = (n- l)(d-i), Q<i<d-1, (6-110) 

where its underlying graph is the cartesian product of d-tuples of cyclic group Z„. Following 
Ref. [1] , the amplitudes of walk and the spectral distribution in the symmetric product of 
graphs can be obtained in terms of sub-graphs. Finally we obtain the following expression for 
the amplitude at origin and spectral distribution 

^ = S n^n(d-0! ^^"~"^ + ^^' 

respectively. 

Also one can show that its idempotents {£'0, Ei, Ea} are symmetric product of d-tuples 
of corresponding idempotents of trivial schemes Kn- That is, we have 

J-, Jn Jn ^ 

-C/o = — <8) — (8) ... ® — , 
n n n 
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E,^ y: {In--)®-®...®- 

permutation 



E,^ Y: (In--)® (In--)...® (In--)®-®...®-. (6-112) 

^ n n n n n 

permutation s / 



Therefore, for the eigenvalues Pjj, and dual ones Qij we get 

PlJ=Ct{C')-\n-lr^K,{^l 



rrii 



Qij = '-^Cf{Ctr\n - ly-'Uj), (6-113) 
where Kk{x) is the Krawtchouk polynomials defined as 

K,ix) ^j^C^C^m-md - (6-114) 

i=l 

and C[ — -j^^fZk)]. Then by using Eq.(2-10) we obtain the amplitude as 

(0,|0o(^)) = # E e-^'''"^^^'^'' Qkr (6-115) 



7 Conclusion 

Continuous-time quantum walk on underlying graphs of the association schemes have been 
studied by using the irreducible modules of Bose-Mesner and Terwilliger algebras connected 
with them, where the irreducible modules of Terwilliger algebra and dual eigenvalues of as- 
sociation schemes play an important role. Also Continuous-time quantum walk on distance 
and strongly regular graphs are investigated by using the spectral distribution associated to 
irreducible representation of associated TerwiUiger algebras. Although the powerful method 
of spectral distribution associated to irreducible representation of Terwilliger algebra seems to 
work for continuous-time quantum walk on distance regular graphs, we expect that one can 
further develop it for some non-distance regular graphs of association schemes and discrete- 
time quantum walk on graphs of association schemes. These are under investigation. Also, it 
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seems that using this formahsm one can study the continuous-time and discrete-time quan- 
tum walk on some graphs which lack association scheme structure, but they have the same 
staratification structure as those of scheme graphs, provided that quantum walk starts from a 
distinguished site (see Ref.[l]). 

Appendix A 

In this appendix, we study the method of symmetrization of group schemes of non-symmetric 
one. If we have a e Cj but is not in Cj, then the association scheme is non-symmetric. 
In order to construct a symmetric group scheme from a give non-symmetric one, we need to 
define the following classes sum 



Ci for i = 0, 

(A-i) 

A + Cr fori = Z + l,...,^, 



where, Ci and Q for i = 0, 1, I are real and i — I + 1, ^^^y^, are complex. 

One can show that the above defined classes sum yield the following relations among 
themselves 



d+l + l 



for i,j = 0, 





\c,\ 


2\G 





v,k=0 AJ^K^) 



d+l+l 
2 



v,r=l+l 



for i = 0, 1, / and j = / -t- 1, 



d+l- 



2 ' 







2\G 
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d+l+l 
2 



(A-ii) 



u,r=l+l Xi/(1) 

for i,j — l+l, ^ii±^. Therefore, the corresponding idempotents {£'0, -Bi, E d+i+i } of group 
association scheme are its irreducible CG-modules projection operators, i.e 



Ek 



Xfe(l) 



d+l+l 



\G\ ^r- 



3 



(A-iii) 



|G| ^3 



EjJ+i(Xfe(a;j) +Xfc(a;j))Cj- 



for A; = 0, 1, I and 



\G\ 



d+l + i 



for A; = Z + 1, 

Obviously the above defined association scheme is symmetric. It is rather straightforward 
to see that its eigenvalues Pj^ and dual ones Qife are 



(A-iv) 



ik 



for A; = 0, 1, / and 



ik 



for A; = 0, 1, 



for i = 0, 1, /, 



^^(XiK)+XiK)) fori = Z + l,...,^ 



^{XiM + XiK)) for z = 0, 1, /, 



(A-v) 



(A-vi) 



Qik 



for A; = 0, 1, / and 



Qik 



dkXk{oti) for i = 0, 1, 

^(XikW + Xik(ai)) fori = Z + l,...,^i|±^, 



for i = 0, 1, /, 



24Xfc("i) 

4(Xit(«i) + Xfc(ai)) for i = / + 1,...,^^, 



(A-vii) 



(A-viii) 



for A; = 0,l,...,^±i±^. 
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In fact, eigenvalues (dual eigenvalues) are sum of real and non-real contributions. 
In section 5, using the above prescription we have studied the continuous-time quantum 
walk on cyclic groups. 

Appendix B 

In this appendix we prove the following lemma in connection with the equality of continuous- 
time quantum walk amplitudes on the vertices belonging to the same associated class. 

Lemma 1. Let qik(t) denote the amplitude of observing the continuous-time quantum 
walk at vertex i e Tk{a) at time t. Then for a class of distance-regular graphs, the amplitude 
Qik is the same for all vertices of associated class k, for all t. 
Proof. 

As the distance-regular graphs form an association scheme, therefore they have corresponding 
TerwiUiger algebra. Also, for Terwilliger algebra T, there exist a unique irreducible T-module 
which has reference point 0, where we called Wq, such that the unit vector Eq.(2-24) is a basis 
for -BfeWo, and the irreducible T-module Wq is orthogonal to its other irreducible T-modules. 
Actually the basis of irreducible T-module Wq can be obtained simply by acting the operator 
A'^ on reference state |0o) repeatedly and also the other irreducible T-modules can be obtained 
by repeaded action of the operator over some set of orthogonal vectors which are orthogonal 
to irreducible T-module Wq and at the same time they vanish under the action of the operator 
A~ . The remaining part of proof is almost similar to the proof presented in Appendix A of 
Ref. [1]. 

Appendix C 
List of some important distance regular graphs 
1. CoIIineeirity graph gen. octagon GO(s, t). 
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^ + 1) + 4(.+i-;£[(l.+v2^) ^(^ - . + 1 + x/2it) + + 1 + 1). 

Intersection numbers: 

Co = s(t +1), Ci = st, (■:2 = .st, C3 = St, 

= 62 = 1, &3 = 1, 64 = ^ + 1- 

2. CoUineairity graph gen. dodecagon GD(s, 1). 

Intersection numbers: 

Co = 2s, Ci = S, C2 = S, C3 = S, C4 = S, C5 = s, 
61 = 1, 62 = 1, &3 = 1, ^4 = 1, &5 = 1, &6 = 2. 

3. M22 graph. 

^^-L-S(x + 4:) + ^d{x + 3) + l^5{x - 1) + \5{x - 4) + 3^^5(0: - 7). 
Intersection numbers: 

Co = 7, ci = 6, C2 = 4, C3 = 4, 

61 = 1, 62 = 1, &3 = 1, ^4 = 6. 

4. Incidence graph pg(A; — l,k — l,k — 1), /c = 4, 5, 7, 8. 

/X = 2F('^(a: -k) + 5{x + A;)) + ^5{x) + - v^) + + v^)). 

Intersection numbers: 
Co = /c, ci = A; - 1, C2 = /c - 1, C3 = 1, 
61 = 1, 62 = 1, by, = k — 1, 64 = A;. 

5. Coset graph doubly truncated biuciry Golay. 

^ = ^^s{x + 11) + §S{x + 3) + i|5(a; - 5) + ^,6{x - 21). 
Intersection numbers: 
Co = 21, ci = 20, C2 = 16, 
61 = 1, 62 = 2, 63 = 12. 

6. Coset graph extended ternciry Golay code. 
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^, = ^^S{x + 12) + + 3) + - 6) + ^,5{x - 24). 

Intersection numbers: 
Co = 24, ci = 22 , 02 = 20, 

= 62 = 2, 63 = 12. 

7. Wells graph. 

^ = ^S{x + 3) + + ^5) + 6{x - ^5) + - 1)) + ^6{x - 5)). 

Intersection numbers: 
Co = 5, ci = 4, C2 = 1, C3 = 1, 
61 = 1, 62 = 1, 63 =4, 64 = 5. 

8. 3-Cover GQ(2, 2). 

// = + 3) + pix + 2) + ^^(a; - 3) + |(5(a; - 1) + j-^S{x - 6)). 
Intersection numbers: 
Co = 6, ci = 4, C2 = 2, C3 = 1, 
61 = 1, 62 = 1, 63=4, 64 = 6. 

9. Double HofFman-Singleton. 

^ = + 2) + 5(x - 2)) + + 3) + 5(x - 3)) + ^,(5(x + 7) + 5(x - 7)). 

Intersection numbers: 

Co = 7, ci = 6, C2 = 6, C3 = 1, C4 = 1, 
61 = 1, 62 = 1, 63 = 6, 64 = 6, 65 = 7. 

10. Foster graph. 

= i(5(x) + l(5{x + 1) + - 1)) + j^(5(x + 2) + - 2)) + ^{5{x + 3) + 5(x - 3)) + 
^(5(a; + ^) + 5(x- v^)). 
Intersection numbers: 

Co = 3, ci = 2, C2 = 2, C3 = 2, C4 = 2, C5 = 1, ce = 1, C7 = 1, 

61 = 1, 62 = 1, 63 = 1, 64 = 1, 65 = 2, 66 = 2, 67 = 2, 63 = 3. 
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In all of the above examples, is defined in terms of bk and Ck as 

CoCi...Cfe_i 

«fe = -TT 7—- (C-l) 

Ol02...0k 
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Figure Captions 

Figure-1: Shows edges through a and /3 in a distance regular graph. 
Figure-2: Shows the graph of symmetric group 5*3. 
Figure-3: The Petersen graph. 
Figure-4: The Johnson graph J(4, 2). 



